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ABSTRACT 

We present ATCA interferometric observations of the old (13 Myr), nearby 
(86pc) classical T Tauri star, PDS 66. Unresolved 3 and 12 mm continuum 
emission is detected towards PDS 66, and upper limits are derived for the 3 and 
6 cm flux densities. The mm-wave data show a spectral slope flatter than that 
expected for ISM-sized dust particles, which is evidence of grain growth. We 
also present HST/NICMOS 1.1 micron PSF-subtracted coronagraphic imaging 
observations of the circumstellar environment of PDS 66. The HST observations 
reveal a bilaterally symmetric circumstellar region of dust scattering ~0.32% of 
the central starlight, declining in surface brightness as r~ 4,53 . The light-scattering 
disk of material is inclined 32±5° from a face-on viewing geometry, and extend to 
a radius of 170 AU. These data are combined with published optical and longer 
wavelength observations to make qualitative comparisons between the median 
Taurus and PDS 66 spectral energy distributions (SEDs). By comparing the 
near-infrared emission to a simple model, we determine that the location of the 
inner disk radius is consistent with the dust sublimation radius (~1400 K at 0.1 
AU). We place constraints on the total disk mass using a flat-disk model and find 
that it is probably too low to form gas giant planets according to current models. 
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Despite the fact that PDS 66 is much older than a typical classical T Tauri star 
(<5 Myr), its physical properties are not much different. 



Introduction 



Classical T Tauri stars (cTTS) are low mass, pre-main sequence stars that are sur- 
rounded by accretion disks, believed to be an integral part of star formation process. Much 
of the cloud materi al from which th e star forms is thought to pass through the disk and ac- 
crete onto the star (Shu et all 1987 ). PDS 66, discussed in detail in this paper, was identified 
as a cTTS by iGregorio-Hetem et al.l (119921 ). The subsequent stage in the pre-main sequence 



evolution scenario (e.g. lLada 1119871 ) is weak-emission T Tauri, when signs of active accretion 



are no longer observed. Infrared excess stud ies in young clusters indicate that disks dissipate 
on ~10 Myr timescales (IHaisch et al.l 120011 ). although there are some examples of disks that 
still contain their primordial gas at or beyond 10 Myr (e.g. TW Hya, CQ Tau, St 34, Hen 
3-600). 

Several processes affect protoplanetary disk evolution, important among them: (1) ini- 
tial conditions in the cloud, such as the reservoir of material, the initial angular momentum, 
and the magnetic field; (2) the co nservation of angular momentum which leads to disk 
spreading as accretion takes place (IHartmann et al.l 119981 ). therefore the size of the disk 
likely depends on its age; (3) grain growth and fragmentation affects the dust opacity, and 
therefore the observed SED; (4) the presence of a planet is expected to create gaps in the 
disk, leading to minima in the SED which cannot be ac counted for by other me ans; (5) disk 



dissipation via stellar winds or photoevaporation (e.g iHollenbach et al.l 120001 ) will dictate 
how long a disk is observable. It is of current interest to study disks over a range of ages to 
test current evolutionary models. 

The combination of infrared and sub-mm observations provide a wealth of information 
about protoplanetary disks. The near-infrared emission above stellar ph otospheric levels, o r 
lack thereof, provides information about the inner edge of the disk (e.g. ICalvet et al.l 120021 ). 
The spectral shape of the silicate feature at 10 /im is a probe of gra in growth and crys 



talline fraction of silic ates that must have formed within the disk (e.g. iKemper et al.l 12004 



Bouwman et al.ll2008l ). The compositional differentiation of the dust population within the 
disk is a relic of the temperature history during epochs of cloud collapse and early disk evo- 
lution. The slope of the far-infrared emission is affected by the amount of disk flaring (radial 
dependence of disk scale height above the midplane). If the disk is geo metrically flatter, then 
the longer-to-shorter wavelength flux density ratio (e.g. 70/im/30/im, Bouwman et al. 20081 ) 
will be relatively lower than if the disk is more flared. Finally, the slope of the millimeter 
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emission is sensitive to grain growth, in that as grains become larger, they will contribute 
more opacity at longer wavelengths. Therefore, an evolving grain size distribution will result 
in a different slope in the millimeter spectrum as grains grow to millimeter sizes and larger. 
The earliest stages of planet formation require growing micron-sized grains into centimeter- 
sized pebbles, and on up to meter-sized rocks. Observational evidence that this process takes 
place is becoming more com mon, and the prese nce of cm-sized grains has been inferred in a 
few circumstellar disks (e.g. IWilner et al.ll2005l ). 



PDS 66 (MP Mus, Hen 3-892), a KUVe cTTS , is a member of Low er Centaurus Crux 
(LCC), a subgroup of the Sco Cen OB association (IMamajek et al.ll2002l ). exhibiting excess 
infrared emission fro m 3 to 70 ittm that has been i nterpreted as arising fr o m a primordial cir- 



cumst ellar disk (e.g.. lGregorio-Hetem et al.lll992l ; ISilverstone et al.ll2006l ). iHillenbrand et al. 



(120081 ) presented the 3 to 70 //m SED, fit to single-temperature photosphere plus disk black- 
body curves (their Figure 6), and interpreted the emission as arising from a wide range of tem- 
peratures. The age of the LCC association has been determined by placing t he known mem- 



bers o f LCC on an HR diagram and using theoretical evolutionary tracks. IMamajek et al. 



( 120021 ) estimated the distance to PDS 66 (86 pc) and the age for LCC (17± 5 Myr) based on 
three evolutionary models. Recent work by iPriebisch fc Mamajek I ( 120081 ). which includes 
more members of LCC, suggests a spatial dependence of the ages, such that the northern 
region appears to be older and the southern part younger. PDS 66 is in the southern part 
of LCC which appears to have an age of ~13 Myr. The periodic photometric variation of 
PDS 66 at visible wavelengths is a signature of stellar rotation, since it is thou ght that accre- 



tion o nto the star produce s observable hot spots th a t co-rotate with the s tar (IBouvier et al. 



19931 : Edwards et al.lll993l : iKundurthv et al.ll2006t >. iBatalha et all (119981 ) find that PDS 66 
has a rotation period of 5.75 ± 0.03 days. PDS 66 is still actively accretin g at an age that 
is substantially older than the median age for cTTS. iPascucci et al.l (120071 ) derived a mass 
accretion rate for PDS 66 of 5xlO _9 M /yr by modeling the Ha line profile; the presence of 
accreting gas confirms that this is still a primordial disk. 

Here we report new multi-wavelength observations combined with information extant 
in the literature to better characterize the disk geometry and properties of the circumstellar 
material in the PDS 66 proto-planetary disk system. This work offers constraints against 
which future models of the system can be tested, and that will eventually reveal the similari- 
ties and differences between PDS 66 and other cTTS. In §2], we describe the observations and 
calibration of our Australia Compact Telescope Array data, the near-IR imaging observations 
using the Near Infrared Camera and Multi-Object Spectrometer (NICMOS) coronagraph on 
the Hubble Space Telescope (HST), followed by a presentation of the SED composed mostly 
of data that we extract from the literature. £J3] describes the disk parameters that we derive 
from the near-IR scattered light observations. $Q describes the constraints on the inner disk 
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properties, derived by fitting models to the NIR region. We discuss our results for PDS 66 in 
the context of other T Tauri stars in £|5l as well as the potential for planet formation within 
the disk. We finish with a summary of our results in §6J 



Observation and Data Reduction 



2.1. NICMOS Data 



Hubble Space Telescope (HST) imaging observations of PDS 66 were obtained as part of 
a Near Infrared Camera and Multi-Object Spectrometer (NICMOS) coronagraphic imaging 
survey (HST/GO 10527; PI: Hines) of a Spitzer-selected sample of 41 stars exhibiting strong 
infrared excesses. The observational s trategy for all targ et stars in the GO 10527 sample is 
identical and is described in detail by lHines et al. r|2oo7h . 



PDS 66 was observed at two celestial field orientations (spacecraft roll angles) to enable 
differentiatin g between artifacts in the point spread function (PSF) from true circumstellar 
features (e.g. lSchneider et al.ll200ll ). The PDS 66 coronagraphic observations (r = 0.3" image 
plane mask with NICMOS camera 2; image scale = 75.8 mas/pixel) obtained on 13 April 
2006^ utilizing two contiguous HST orbits. In each orbit the target was autonomously posi- 
tioned into the coronagraph following short (0.915 s) unocculted target acquisition images in 
the F171M filter. These images were also used in conjunction with the spacecraft positioning 
information to determine the position occulted (< 5 mas RSX differential uncertainty). Each 
orbit resulted in 2591 s of total coronagraphic integration time with the F110W filter. At 
the end of each orbit, PDS 66 was directly imaged (non-coronagraphically) with four dither 
offset 6 s exposures to enable flux density normalization for PSF subtraction. 

Images d erived of PDS 66 and t emplate PSFs (discussed below) were created following 
the method of ISchneider et al.l (120051 ) (SSH05) and references therein. Within each visit, the 
individual images were median combined after verifying their repeatabilities. A suite of PSF 
subtracted PDS 66 images were created from the target and candidate PSF template median 
combined images from an ensemble of bright GO/10527 spectral type G-F targets. Initial 
target:template flux scaling was performed by using 2MASS catalog J (a close surrogate to 
NICMOS F110W) magnitudes. 



After flux scaling and astrometric registration (see SSH05), all suitable PSF subtracted 



1 Visit numbers 19 and 20 as per the Multimission Archive at Space Telescope 
http: / /archive, stsci.edu/hst/search.php 



- 5 - 



imaged showed similar: (a) scattered light excesses at r < 1.0", (b) oversubtractions at 
r>1.2". The latter is most likely result from chromatic effects in PSF subtraction. With a 
J-H color index of +0.64, PDS 66 is quite red compared to ~ +0.2 for G-F main sequence 
stars used as initial PSF templates. 

Better PSF subtractions were obtained using a high SNR, well color-matched, coro- 
nagraphic calibration PSF template derived from observations of the calibration star HD 
82558 (spectral type K0, J = 6.08, J-H = +0.48) observed in HST GO program 10177 (PI: 
G. Schneider). The two PDS 66 minus HD 82558 PSF subtracted images were median com- 
bined after rotating difference images to a common celestial frame and masking in each a 
small local artifact as well as the HST diffraction spikes. The scattered-light results discussed 
in this paper are derived exclusively from the two-orientation recombination of these higher 
fidelity PSF subtractions. 



2.2. ATCA Data 

The Australia Telescope Compact Array (ATCA) was used to observe the continuum 
emission towards PDS 66 in the 3 mm (94.1 GHz), 12 mm (18.9 GHz), 3 cm (8.6 GHz), 
and 6 cm (4.8 GHz) bands between 28 August and 2 September 2005. The phase center 
was 13 h 22 m 8.53 s in RA and -69° 38' 12.17" in DEC (J2000), which is offset by 5.1" in RA 
from PDS 66 to avoid artifacts created at phase center. The ATCA interferometer consists 
of six 22-meter telescopes. We used the H214 configuration, where 5 antennas are arranged 
with baselines between 82 and 247 m, and the sixth antenna is on a 4.5 km baseline. The 
atmospheric phase stability was too poor to calibrate the longest baseline, and PDS 66 is 
too faint for self-calibration, and therefore data from the sixth antenna were discarded. The 
system temperature for the single sideband receiver ranged from 380-500 K at 3mm, and 
30-60 K at 12mm, 3cm, and 6cm. The correlator was configure d for 32 channels across 128 



MHz of bandwidth in each band. We used the MIRIAD package (ISault et al.lll995l ) to reduce 
and calibrate our data, and all of the calibration sources are from the 1997 ATCA Calibrator 
Source Catalogue by J. E. Reynold^. 

The bandpass, flux, and gain calibrators for the observations in each of the four bands 
are listed in Table [1] with other key observational parameters. Observations in the 3 mm 
band were obtained on four separate days. The map presented in Figured] was produced from 



2 Except those using HD 61005 possessing a bright light-scattering disk, and HD 25427 rejected due to 
instrumental image degradation. 

3 ftp://ftp. atnf.csiro.au/pub/atnfdocs/guides/at. cat 
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one day of the observations which had the best atmospheric phase stability. We measured 
a flux density of 2.2 Jy at both frequency windows (93.1 and 95 GHz) for 1057-797 on 2 
September 2005 using Uranus as the flux calibrator. We assumed 1057-797 had a constant 
3 mm flux over the four days. It was used as a secondary flux calibrator for the other three 
days. The flux scale in the 12 mm band was set by averaging the flux of 1352-63 over two 
runs, which lowered our uncertainties. 



Table 1: ATCA Observing Parameters 





3 mm 


12 mm 


3 cm 


6 cm 


Primary Flux Calibrator 


Uranus a 


1934-638 


1934-638 


1934-638 


Assumed Flux Density (Jy) 


4.5 


0.99 


2.8 


5.8 


Phase Calibrator 


1057-797 


1352-63 


1251-713 


1251-713 


Derived Flux Density (Jy) b 


2.2±0.3 


0.98±0.03 


0.94±0.03 


1.02±0.03 


Bandpass Calibrator 


1921-293 


1253-055 


1253-055 


1253-055 


PDS 66 










Derived Flux Density (mJy) 


22.0±3.3 


0.41±0.06 


<0.20 d 


<0.25 d 


Internal Uncertainty (mJy) 


2.4 


0.06 






Calibration Uncertainty (mJy) 


2.3 


0.01 






UV Range (kA) 


21-69 


3-16 


1.5-7 


0.5-4 


FWHM Primary Beam (') 


0.5 


2.5 


5.5 


9.9 


FWHM Synthesized Beam (") 


3.1x2.2 


10.2x9.2 


23.8x20.2 


43.5x34.9 



a The flux of Uranus was derived based on an assumed brightness temperature and a uniform disk. 
The uncertainties were estimated by measuring the flux at different elevations during the track. 
c The bandpass is the frequency-dependent gains 
d 3a upper limits 

Table [1] lists the derived flux densities or the upper limits for PDS 66 in each of the 
four bands, along with the uncertainties (internal, calibrated, and total). The signal-to-noise 
ratios for 3 and 12 mm are 6.7 and 6.8, respectively. There was no detection at 3 or 6 cm, and 
so we list their 3-cr upper limits. The fluxes were measured in the UV-domain using UVFIT 
with a point source model. Images were made using INVERT with a robust weighting of 2 
to optimize sensitivity, and the CLEAN algorithm was used to deconvolve the images. The 
restored maps in Figure 1 are shown for 3mm (left) and 12mm (right), with the fwhm of 
the beam shown in the bottom left corner. The 3mm data provide the maximum angular 
resolution of 2.2" or 190 AU at the distance of PDS 66 (86 pc). The visibility amplitudes as 
a function of projected baseline length is shown in Figure [2j 
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Fig. I. — 3mm (left) and 12mm (right) emission from PDS 66. The fwhm of the synthesized 
beam is shown at the bottom of each plot in solid black. The contours are drawn as no", 
where n=-6,-5,-4,-3,3,4,5,6 and <7=3.2 mJy/beam for 3mm and 0.06 mJy/beam for 12mm. 
The negative contours are dashed and the positive contours are solid. 



- 8- 




20 40 60 

UV distance (kX) 




5 10 15 

UV distance (kX) 



Fig. 2. — Visibility amplitude as a function of projected baseline length for PDS 66 observa- 
tions in the 3mm (left) and 12mm (right) bands. The dotted line is the expected amplitude 
for zero signal, and is calculated in MIRIAD based on the observed variance within each 
bin. The data points in each bin are the azimuthally averaged visibilities. Since the profile is 
consistent with being flat, the disk of PDS 66 is not spatially resolved in these observations. 
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2.3. Spectral Energy Distribution 



The SED for PDS 66 shown in Figure [3] was derived from data obtained with several 
in struments, taken at d ifferent times. The U,B,V,R C ,I C photometry are taken from the mean 
of iBatalha et al.l (119981 ). since they show that PDS 66 varies by 0.1 mag over 5.75 days at 
these wavelengths. We include this photometric variability as a source of uncertainty in 
the optical photomet ry, and find that the fluxes agree with the B and V p hotometry from 



the T ycho-2 catalog flHog et al.ll2000bh . We use 2 MASS J,H,K photometry (ISkrutskie et al. 



20061 ). and IRAS 12 and 25 /im photometry. We did not include the IRAS 60 /im data 
point because the field was contaminated by cirrus. The Spitzer data, including IRAC, 
MIPS, and a low-resolution IRS spectrum, wer e obtained and reduced by the Formation and 
Evolution of Planetary Systems (FEPS) team (IMeyer et al.ll2006l ; iBouwman et al.ll2008l . and 
references therein) . Th e 1.2 mm point comes from the Swedish-ESO Submillimeter Telescope 
( [Carpenter et al.ll2008l ) and the 3 and 12 mm data are those described in §2.21 



We estimate the extinction towards PDS 66 by comparing t he der eddened spectrum to a 
standard spectrum for a K1IV star from lKenyon fc Hartmannl (119951 ) normalized to I-band. 
We found that extinction corrections derived from color excesses in (B-V) overestimates the 
fluxes at A longer than I-band, and underestimates the fluxes at A shorter than I-band; 
while color excesses in (J-H) shows the opposite trend. We find the best agreement for 
extinction as estimated from color excess with (R-I) c and (V-J), which are in good agreement 
with each other . We adopt v4„=0.7±0.2 mag. This is different than the value derived in 
Mamajek et al.l (120021 ). where they used the weighted mean of (B-V) and (V-J) to derive 
v4„=0.17±0.07. The reason for the discrepancy might be due to the slightly elevated B- 
band fluxes produced by accretion onto the star, instead of arising from stellar emission 
alone . We update the lumi n osity accordingly to 1.1 L Q , utilizing the bolometric correction 
from iKenyon fc Hartmannl (119951 ) for a K1V star. This h igher luminosity decreases th e 
age of PDS 66, assuming the same temper ature (5035 K ) as IPriebisch fc Mamajek I (120081 ). 
Using the stellar evolutionary models of iBaraffe et al.l (119981 ). we estimate an age of 11 
Myr. These same mode ls were used to estimate the age of the moving group (13 Myr) to 



which PDS 66 belongs (IPriebisch fc Mamajek 1 120081 ). this mean group age is more robust 



than the age derived from the position of PDS 66 on the HR-diagram alone, the results are 
consistent within t he uncertainties . We corrected the flux densities for extinction using the 
extinction law from iMathis I (Il990l ) translated into the Johnson- Cousins color system. Color 
corrections were applied to the IRAC (IRAC Data Handbook, v. 3.0) and MIPS (MIPS Data 
Handbook, v. 3.2) photometry. Color corrections were not applied to the IRAS data because 
the correction factors were close to unity (within 3%). The total uncertainties represent the 
internal and absolute calibration errors, and extinction corrections where appropriate. The 
literature data are listed in Table [21 and the dereddened SED is shown in Figure [3J 
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The SED of PDS 66 is compared to the median Taurus SED (ID'Alessio et al.lll999l ) in 
Figure [3], where we have normalized to the H-band. One should recall that: (1) the SpT 
of the Taurus sources used to construct the median are K5-M2, whereas the SpT of PDS 
66 is Kl, and (2) PDS 66 is likely more face-on than the mean expected inclination of the 
Taurus median (which is probably closer to 60° assuming a random distribution of inclination 
angles). With these caveats, we note the flux deficit relative to the Taurus median in the 
4-20 /im region. This could be the result of disk-clearing close to the star, or shadowing of 
parts of the disk further out in radii by the puffed-up rim. A detailed analysis accounting 
for the inclination and SpT is described in §4.11 Another feature which stands out from the 
Taurus median is the flux deficit near 70 /im. Since PDS 66 is more face-on than the expected 
in clination of the median Taurus disk, a higher FIR luminosity is expected (see Figure 5A 



in 



Dullemond et al.ll200ll ). On the other hand, the SpT is earlier for PDS 66, resulting in 
a larger scale height of the PDS 66 rim, so more s hadowing of neighborin g parts of the 
disk would lower the FIR emission (see Figure 5C in iDullemond et al.l 120011 ). However, the 
difference between the PDS 66 and the Taurus median SED may be a real physical signature 
of a less flared disk with larger grains. Discussion of these possibilities is presented in §5.11 



A(/xm) F v (Jy) o totai (Jy) 



0.37 


0.031 


0.003 


0.44 


0.122 


0.006 


0.54 


0.263 


0.007 


0.64 


0.392 


0.011 


0.80 


0.546 


0.010 


1.24 


0.791 


0.023 


1.66 


0.909 


0.018 


2.16 


0.789 


0.013 


3.54 


0.657 


0.014 


4.51 


0.521 


0.012 


7.74 


0.471 


0.010 


12.0 


0.882 


0.053 


23.7 


1.874 


0.120 


25.0 


2.071 


0.145 


71.5 


1.54 


0.12 


156.0 


2.06 


0.29 


1200 


0.207 


0.02 



Table 2: Reddened fluxes for PDS 66 from the literature, as described in the text. 
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Fig. 3. — Dereddened spectral energy distributi on of PDS 66 compose d of data from the 
literature, the low-resolution IRS spectrum from iBouwman et al.l (120081 ). and new millime- 
ter points from ATCA des cribed in §2.21 A K1IV star with T e ^=5035 K (colors from 
Kenyon fc Hartmann ( 19951 )) and normalized to I-band is shown as the dotted line. The 



Taurus mean SED normalized to H-band is sho wn for reference as the solid black line, with 
the upper and lower quartiles in dash-dot lines (ID'Alessio et al.lll999l ). 
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3. Analysis and Results From the Scattered-Light Disk 

The HST/NICMOS 1.1 fim scattered-light image of the PDS 66 circumstellar disk is 
shown in Figure H] (right image). A faint, point-like, object apparently superimposed upon 
the disk at a celestial PA of 359°, 1.2" from PDS 66 (F110W mag of ~17.2; 5-F110W ~+9.0) 
appears in the NICMOS image. VLT/NACO J, H, K images obtained in May 2003 also reveal 
the object, but not the disk (M. Janson, priv. comm. May 2006). The VLT images were 
used in combination with the NICMOS observations and establish through differential proper 
motions measures that this is a background star and not a faint companion associated with 
PDS 66. 

The morphology of the 1.1 /zm PSF-subtracted two-orientation combined scattered-light 
coronagraphic image of the PDS 66 disk suggests an inclined disk with a PA=169±6° for 
the disk major axis. The disk is robustly detected on the major axis to r~1.5", with a low 
surface brightness scattering component extending to r<2" (~170 AU). The disk appears 
mirror symmetric (at low spatial frequencies) about the minor axis. The eastern hemisphere 
of the disk (to the left in Figure H]) is significantly brighter (at equal radii) then the opposing 
side of the disk at all diametrically opposed azimuth angles, suggesting preferential for- 
ward scattering toward the morphological minor axis at PA~79°. The azimuthal brightness 
asymmetry with respect to the apparent major axis of the disk is indicative of an inclined 
disk. Assuming the apparent eccentricity is due to projection, the inclination is estimated 
to be 32±5° from face-on based upon the apparent major and minor axis lengths (mean 
axial ratio=l:0.848) determined from five equally spaced isophotal fits from r=12 to 22 pixel 
(0.91<r<1.67"). 

We measure the 1.1 fim scattered light flux density of the PDS 66 circumstellar disk 
from an effective inner working angle of 0.41" (one resolution element beyond the 0.3" corona- 
graphic mask radius) to 3.5" (beyond which no significant disk-scattered starlight contributes 
to the total brightness) as 2.7±0.4 mJy This excludes a small area circumscribing the coron- 
agraphic obscuration (shown as black in FigHJand rejected as unreliable due to imperfections 
in PSF subtraction very close to the star). Thus, the total disk flux density may be some- 
what greater. We adopt an F110W absolute instrumental sensitivity of 1.26 xl0~ 6 Jy per 
ADU/sec/pixel and a Vega system photometric zero point of 1775 Jy for a spectrally flat 
source. By scaling the brightness of the HD 2558 template to the PDS 66 observations, we 
find the F110W magnitude of PDS 66 to be ~8.29 (= 860 mJy). Thus, the fraction of 1.1 
/im starlight scattered by the disk at r>0.41" is approximately 0.32%. 

While the PDS 66 disk appears bilaterally symmetric (implying a disk inclined to the line 
of sight), we measured, and fit, the disk's radial surface brightness profile medianed over all 
azimuth angles (FigureHJ). We find that within 0.8<r<1.5", the disk's azimuthally medianed 
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Fig. 4. — NICMOS F110W 2-roll combined PSF-subtracted imaging observations of the 
PDS 66 circumstellar disk. Left: Azimuthally medianed radial surface brightness profile 
(excluding unreliable data close to the star, shown masked as black in plot). Error bars 
are standard deviations (la) about the median surface brightness (SB) at corresponding 
radii. Dotted line is extrapolation of a power law fitted in the region in blue (0.8" <r<1.5") 
extrapolated to r=3.5". Plotted below is the SNR corresponding to radial flux density 
measurements in each one pixel wide annular zone. Right: 7.5x7.5" field (north up, east left) 
centered in the NICMOS coronagraph; loglO stretch from +0.001 to +10 ADU s _1 pixel -1 (1 
ADU s" 1 pixel -1 = 218 /xJy arcsec -2 ). Filled gray circle (r=0.41") is one resolution element 
beyond image plane obscuration. 
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3.1. Grain Properties Inferred from Millimeter Wavelengths 



The growth of grains may be inferred from the slope of the millimeter spectrum, where 
the emission goes as u a . Assuming the Ray leigh- Jeans limit, and that the opacity is a simple 
power-law going as it follows that a = 2 + (3 if the emission is optically thin. Smaller 
grains (sizes <C A, e.g. ISM-like grains) have /3 ~ 2, while larger grains (sizes ^> A) have (3 ~ 
and the emission follows the shape of a blackbody. Potentially complicating this simple 
SED analysis, cTTS occasionally show signs of ionized stella r winds, which can dominate the 
centimeter e mission with spectral s l opes ranging from if~ 0A (iMezger et al.lll967l) for optically 
thin to z/ 6 flPanaria & Felli Ill975l : Iwright & Barlow Ill975l : blnon Ill975h for optically thick 
emission. The stellar wind can contribute significantly to emission in the millimeter, so 
corrections for a wind component must be applie d so that only therma l dust emission is 
considered in deriving limits on the spectral slope. iRodmann et al. fit the centimeter 



emission for young stars in Taurus, and found a range of slopes from v~ 0A to z/ a4 . We use 
our 3a upper limits fro m the centimeter observations with the adopted range of slopes from 
Rodmann et al.l (120061 ) to estimate possible contributions from stellar winds. 



Figure [5] shows the mm- wave spectrum, with a linear least-squares fit to the data points 
(dotted black line), and estimates of the stellar wind (dashed and dotted lines) going through 
the upper limits to the centimeter data. Although we do not find evidence for a strong stellar 
wind, we point out that the values extrapolated from the 3a upper limits could contribute 
significantly at 12mm (between 43-65%). However, since we did not detect such emission at 
3 or 6 cm, we assume that stellar wind is negligible. 

We measure a slope of 2.4±0.1 for the millimeter emission, which corresponds to a j3 
of 0.4±0.1. This derivation of f3 assumes purely optically thin emission, yet there can be 
a non-negligible contribution from optically thick emission which would cause (3 to appear 
lower. The correction to J3 to account fo r some optically thick emission, defined using A, 
is described in detail in iBeckwith et al.l (119901 ) (see their equation 20). We follow their 
same definitions for the temperature, surface density, and opacity profiles; and their same 
normalization of 0.1 cm 2 /g at 10 12 Hz (0.3 mm) for the total (dust+gas) opacity. 

A was derived assuming that the density and temperature decrease with radius as r~ p 
and r~ q , respectively, and so it depends on the density and temperature profile indices, p and 
q. A is we akly dependent t o R^. fc, M disk , and i as they are captured in the lnf factor (Equa- 
tion 16 in IBeckwith et al.l (1990)). Since we only have crude estimates of these quantities, 
we can only approximate the contribution to (3 from optically thick emission. We gener- 
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ousl y estimate A assuming fti,2mm0, P=l-5, g=0.75, i?di S &=170 AU@ (see $3]), M^ s fc=0.005 
M Q ([Carpenter et al.ll2005l ). and i = 32° (see jj3]). The v alue for p is at the high end of a 
reasonable range according to lAndrews fc W illia ms I (120071) ; and q is from the optically thin, 



geometrically flat di sk discussed in 



lower than this (e.g. iBeckwith et al. 



Adams et al.l (119871 ). although typical values for q are 



1990h . We find A -0.2 and (3 COTr < 0.5. 



This low value of 13 may indica te particle sizes near the mm-regime according to the mod 
els of iMiyake fc Nakagawa I (119931 ). where they show that (3 <1 for a grain size distribution 
of n(a ) oca -3 5 and a maximum grain size of 1mm. More recent models by iD'Alessio et al. 



( 120011 ) also suggest that such low values of (3 are found if the maximum grain radii are in the 
millimeter to centimeter regime. Since the millimeter emission traces the bulk of the disk 



mass, it is possible to infer the disk mass from the flux at these wavelengths. ICarpenter et al. 



( 120051 ) derive a dust mass for the disk of PDS 66 of 5 x 10 5 M , and our 3 and 12 mm fluxes 



give very similar masses even with our new estimate of (3. Note that the derived total dust 
mass may be somewhat higher from the low surface brightness diffuse scattering seen beyond 
75 AU in scattered light. 



4 Ki.2 mm maximizes the correction, since K3 mm or Kvimra lead to a lower value for A. 
5 Higher Rdisk values will lead to lower values of A. 
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Fig. 5. — Long- wavelength spectrum showing the fit to the observed data points (dotted 
line). We measure a slope of -2.4±0.1, which corresponds to a of 0.4±0 . 1. Th e solid and 
dashed lines represent typical slopes for the stellar wind ( Rodmann et al. 20061 ) that has a 
frequency dependence u m , and are constrained by our 3 and 6 cm upper limits (arrows). We 
assume that the effects of stellar wind are probably insignificant. (3 has been corrected for 
optically thick emission, using values for the temperature and density power-laws such that 
the correction is maximized. We find f3 COTr < 0.5. 
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4. Analysis and Results From SED Models 

Detailed modeling of the SED combined with resolved imagery at one or more wave- 
lengths can lead to a better understan ding of the physical structure of the circumstellar 



environment (e.g. ISchneider et all 120031 ) . In theory, one can learn about the composition 
and size of the grains, the disk geometry, and mass surface density as a function of radius. 
However, we are limited in the parameters we can derive since the disk was not resolved in 
the millimeter. For example, the outer disk radius, surface density profile, and disk mass 
are degenerate parameters, since we cannot distinguish between a steep mass distribution 
concentrated tightly around the star and a shallow mass distribution extending further out. 
However, we can still constrain specific physical parameters in regions of the SED which 
are not strongly dependent on degenerate model parameters. We focus on estimating the 
inner disk radius, inferring the degree of flaring, and constraining of the surface density and 
disk mass. We estimate the inner dust truncation radius by fitting a single temperature 
blackbody to the near infrared (NIR) emission (2-5 /mi), described in §4.11 The fixed model 
parameters throughout this analysis are the stellar mass (1.3 M Q ), effective stellar temper- 
ature (5035 K), and the distance to the system (86 pc). We choose between models in §4.11 
using the x 2 statistic. Penultimately, we use our millimeter observations to place constraints 
on the surface density and disk mass in §4.21 using a flat disk model. 



4.1. Inner Disk Radius 

A large fraction of the flux in the NIR can be dominated by emission from the inner 
rim (significantly inside the inner working radius of the NICMOS coronagraph which is 25 
AU at 86 pc), a region where the local disk opacity changes dramatically. The change in 
opacity is caused by the boundary between temperatures that are above and below dust 
sublimation, and the l ack of dust inside of this boundary results in a drop to zero opacity 



(IMuzerolle et al.ll2003l ). The inner edge of the dust disk is in full view of the star, and 
therefore, the temperature is higher than it otherwise would have been if the dust disk were 
not truncated. The gas in the disk near the point of dust sublimation is also hotter, and so 
the disk scale height increases there because the pressure is greater. But, the disk might also 
appear to have a region cleared of dust if the grains have grown more in the inner regions, 
since larger grains have lower opacities. Faster growth is expected closer in because the 
particles have higher Keplerian velocities and the surface density in the disk is higher, thus 
higher collision rates are expected there. A planet forming in the disk would also decrease 
the opacity in its vicinity because it would sweep up the surrounding material, creating a 
gap. 



-18- 



It is interesting to test if any of the above scenarios are consistent with the observations 
of PDS 66, which is possible by estimating the temperature (and thus location) of t he inner 
rim. We calculate its emission as F u =B u xa,rea. (following iDullemond et all 120011 ). where 
the area is a function of the inclination (greater surface area of the inner rim wall is seen 
for more edge-on disks) and vertical s cale height of the rim . The U, B, V, R, I, J, H, 
and K stellar colors are estimated from lKenyon fc Hartmannl (119951 ) assuming a K1V star; 
and we extrapolate the fluxes beyond K-band assuming F u oc u 2 . The central luminosity, 
inclination, and vertical scale height are fixed at 1.1 L Q (see §2.31) . 32° (see §3]), and 0.1 
Rrimi respectively. Since a constraint on the vertical scale height is not available, it is the 
dominant source of uncertainty i n R T i m . Varying the vertical scale height from 0.1R rim to 
0.25R r i m (IDullemond et al.ll200ll ) leads to a change in R r im by a factor of 1.5. 



We fit the temperature and solid angle of the inner rim for the 1-5 /im photometry only, 
assuming the rim emission is a Planck function. We do not expect the 8 fim point to fit well 
with a single temperature blackbody at the inner rim, as it is likely contaminated by silicate 
emission and cooler dust temperatures. The left plot of Figure [6] shows the fit for the rim 
temperature and radius which give the minimum xt'- T r i m =1277 K and R r i m =0.07 AU. The 
right plot of Figure O shows the reduced x 2 as a function of the rim temperature, where the 
temperature grid explored exceeds the range shown on this plot (900-1400 K). We find that 
the rim temperature is within 1190-1370 K for a 3a confidence level, and the probability of 
these data being drawn from a model with this temperature is 24.5%. This result indicates 
that the dust temperature and inner rim radius are consistent with being located at the dust 
sublimation radius. 
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A(yu,m) Rim Temperature (K) 



Fig. 6. — The left plot shows the best-fitting rim temperature (1277 K) and radius (0.07 
AU) for a single temperature blackbody (solid line) model. The black diamonds are the 
data points for the excess (disk) emission: J, H, K, 3.6, and 4.5 jum, and their error bars 
include the internal, calibration, and extinction uncertainties. Note that we did not include 
the 8 (Ltm data point in the fit. Plotted on the right is the reduced x 2 as a function of rim 
temperature, showing that the temperatures within 3a that are fit by the model. 
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4.2. Flat Disk Model 



In the millimeter regime, it is a reasonable approximation to model the disk as geomet- 
rically flat because, at these wavelengths, the emission is optically thin and we can observe 
the midplane of the disk where the particle density is the highest. We are interested in 
learning about the possibility of planets forming in this disk, and particularly if the surface 
density and disk mass are sufficient. The following equations can be used to estimate the 
surface density in the planet forming region, where the normalization is at 5 AU. 

We assume a radial temperature distribution with power-law index, q. We take the 
temperature from the best-fitting single-temperature blackbody to the NIR emission de- 
scribed in §4.11 T rim =1277 K at 0.07 AU, and extrapolate to larger radii assuming q=0.5. 
This value of q was as sumed by fitting the SED between 11-70 /im following Equation 10 in 
Beckwith et al.l (Il990l ). The total surface density distribution is given by 



Sfr) 



5AU 



(1) 



We assume p = 1 for the bulk of this discussion, but consider values for p ranging between 
and 1. The assumption of p = 1 and the range in p is motivated by millimeter observations 
that suggest the mean value of p obse rved in disks is close to un i ty and most disks fa ll 
within the range 0<p<l and the (e.g Andrews fc Williams 1 120071 ; iKitamura et al.l 120021 ). 
The surface density normali zation, Sr, is e valuated at 5 AU; and for reference, the minimum 
mass solar nebula (MMSN, lHayashilll98ll ) is £5=150 g/cm 2 for gas and dust. Although the 
value for the MMSN assumes p = 3/2, which is large compared to observed values of p, the 



normalization determined by lHayashil (Il98ll ) is substantially larger than is constrained by 
our observations. 

The quantities described above can be used to estimate specific parameters related to 
the emission coming from a disk. In particular, the flux for a geometrically flat disk is given 
by 



d 2 .1 r 



Rout 



BJ1 



)27rrdr 



(2) 



B v is the Planck function, r v is the optical depth, /i=cosi, and d is the distance to the source. 
Since we are considering emission at wavelengths beyond 1 mm, we can assume the Rayleigh- 
Jeans limit (T <30 K), and that the disk is optically t hin: (1 - e~ Tv ^) ~ t v /j i=k u 'E(R)/ ijl. 



The opacity, k u , is the same equation as described in iBeckwith et al.l (Il990l ). We use an 
inner disk radius of 0.07 AU from the fit to the NIR. With the given assumptions, we can 
solve for the surface density at 5 AU: 

F y c 2 d 2 (2-p-q) 



(3) 
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Figure [7] shows S 5 as a function of outer disk radius, ranging from 1 to 200 AU, calculated 
with Equation [31 and using the 3mm flux. We consider a range for the power-law on the 
surface density between and 1. The vertical line labeled 'r„j S > 10 4 yrs' is the size of a 



disk which has been growing due to viscous evolution for 10 4 yrs (jHollenbach et al.ll2000i ). If 
material is moving inwards, and eventually accreting onto the star, the disk spreads out on 
a viscous timescale (r V i s = r 2 /a) to conserve angular momentum. The region to the right of 
the line represents disks that have higher viscous timescales, and therefore larger outer disk 
radii. We choose a viscous timescale that is much shorter than the age of the star because 
there are other processes that may be interfering with disk spreading, making the viscous 
timescale longer. Thus, we are setting a conservative lower limit (10 AU) to the outer disk 
radius, but we realize that the disk must be larger than this based on both the age of the 
star and the HST/NICMOS observations (§3]). The line labeled 'mm-beam' is the maximum 
beam size from the 3mm observations (130 AU), and we assume the disk size must be smaller 
since the emission from PDS 66 is unresolved ( §2.21) . It is clear in Figure [7] that as the outer 
disk radius marches down, more mass must be packed into a smaller and smaller disk, and 
so the surface density at 5 AU increases rapidly. If we consider the case of p = 1 (dash-dot 
line) we require that 3 < S 5 < 11 g/cm 2 based on the radius constraints. If the disk follows 
a shallower density profile, shown with the dashed and dotted lines in Figure [7J, the spread 
between upper and lower limits on £5 increases. 

The disk mass is tied to the surface density via the equation M = j disk 27rr12(R)dr, and 
so a higher surface density leads to a higher disk mass. If we integrate over the disk, the 
total disk mass is given by 

^ _ 2nE 5 rURl: t P - R^ P ) (4) 
2 — p 

The total (dust+gas) disk mass derived by ass uming p = 1 is . 0004 < M d < 0.001 M . 



For comparison, the total disk mass derived in ICarpenter et al.l (120051 ) from their 1.2 mm 
observation is 0.005 M & . This is approximately 17 M e in dust mass, and 5.3 Mj Up it er in 
total mass (assuming a gas-to-dust ratio of 100). 

Since it is more realistic that the power-law on the temperature and density profiles 
vary with disk radius, we describe how our results change if these indices are different. One 
extreme may be that the disk is optically thick and geometrically thin (q = 0.75) for all 
radii, making the derived disk mass higher by a factor of 2-5. If instead, we keep the same 
temperature profile (q = 0.5) but consider flatter density profiles (by changing the S r profile 
from p = 1 to p = 0), then the limits on the disk mass will increase by a factor of 1.25-2 for 
the range in radii we explore. 
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Fig. 7. — Total (dust + gas) surface density at 5 AU as a function of outer disk radius, 
calculated according to Eqn [3X using our 3mm flux. Shown are the results for three values of 
the power law on the surface density profile: p = 0, p = 0.5, and p = 1. The vertical lines are 
the maximum disk size assumed from the mm observations and the lower limit on the outer 
disk radius is shown as the vertical line labeled l r vis > 10 4 yrs', based on the size of the disk 
estimated from viscous evolution. Assuming the p = 1 model is correct, which is reasonable 
for a T Tauri star, we estimate the surface density at 5 AU and use it to constrain the total 
disk mass to 3< S 5 <11 g/cm 2 and 0.0004 < M d <0.001 M Q , respectively. 
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Discussion 



5.1. How Unusual is PDS 66? 



Despite its age, PDS 66 (~13 Myr) is analogous to the well-stu died cTTS in Taurus (1-3 
Myr). Compared to the sizable sample of stars in Taurus studied by lFurlan et al.l (120061 ) . the 
spectral indices computed between 13 and 25 fim vs. 6 and 13 /zm for PDS 66 (ne-i3 = —0.13, 
^13-25 = 0.36) are typical of the flared, Class II sources (see their Figure 11). Figure [3] also 
illustrates the similarity between the SED of PDS 66 and the median SED of T Tauri stars 
in Taurus. Most sources older than 10 Myr do not display primordial signatures, yet PDS 66 
appears similar to younger cTTS. There may be hints that PDS 66 is slightly more evolved 
than the mean disk in Taurus, but proper accounting for biases between the two SEDs 
must be taken before drawing conclusions. We used a simple model fit to the NIR emission 
discussed in §4.11 to show that the inn er rim is similar in temperature and locat ion to values 
for cTTS derived in other works (e.g. iMuzerolle et al.ll2003l ; iMeyer et al.l 119971 ) . The outer 
disk radius assumed from the HST/NICMOS scattered light i mage (~170 AU) is consisten t 



with the mean outer disk radius derived from mm-observations ([Andrews &: Williams 1120071 ) . 
Although, mm- and NIR scattered light-derived outer disk radii tend to diff er. The outer disk 
radius of IM Lup, for example, is 400 AU derived fr om NIR scattered ligh t (IPinte et al.ll2008l ) 
and only 105 AU from resolved mm-observations (ILommen et al.l 120071 ) . This discrepancy 
is probably due to different sensitivities to mass surface density and emission that traces 
different grain populations. 

We now compare PDS 66 to other sources of comparable in age. TW Hydra, a 10 
Myr old cTTS, appears to have an inner hole that is 4 AU from the star, well beyond the 
dust sublimation radius, inside of which, there is only optically thin material [r v < 0.05, 
Calvet et al.ll2002l). This inner h ole might be created by a planet sweeping up material from 



the disk. ISetiawan et al.l (120081 ) claim that radial velocity variations in TW Hydra could 
be caused by a planet, but subsequent observations and ana lysis suggest that the variations 
are instead c aused by co o l star spots (IHuelamo et al.l 120081 ) . The disk mass of TW Hydra 
(~0.009 M , iKrist et al.l (boooh ) is within a factor o f 2 of PDS 66 (0 005 M Q ), however the 
mass accretion rate of TW Hydra (5 x 10~ 10 M^/yr. IMuzerolle et al.l (120001 )) is lower by an 
order of magnitude. The mass accretion rate may be a means o f discriminating between 
cTTS and transitional disks, as is discussed in lNajita et al.l (120071 ). where transitional disks 
have lower accretion rates for a given mass. TW Hydra falls in the transition disk territory, 
having a lower mass accretion rate for its mass. Interestingly, PDS 66 falls on the border in 
between cTTS and transitional disks for single stars, hinting that it may be evolving towards 
the transition disk population. Perhaps the differences between PDS 66 and TW Hydra are 
related to the initial conditions, such as available circumstellar material. That could explain 
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why TW Hydra appears to be more evolved at a similar age, at least in terms of the inner 
disk properties and accretion rate. 



Andrews fc Williams I (120071 ) find that the mean surface density at 5 AU for their sample 
disks in Taurus (ages <7Myr) is ~14 g cm" 2 , which is comparable to the values we estimate 
here with p = 1 and q = 0.5. The total disk mass constrained by our millimeter observations 
is on the low end of the distribution for disks in Taurus, which may be a sign of disk 
dissipation over time. The value we derive for /3 is smaller than the mean for Taurus, where 
(3 ~ 1, and therefore PDS 66 may have a larger grain population, on average. TW Hydra 



Wilner et al. 


2005; 


Testi et al. 


2003) 



Bouwman et al.l (120081 ) give a detailed analysis of the low-resolution IRS spectrum for 
PDS 66, and they find evidence for grain sizes up to 5 /im. Their observations are probing the 
surface layers of the disk, where smaller grains spend most of their time. Evidence of even 
larger grains are probed by our millimeter observations, but these gr ains are found primarily 
in the midplane of the disk. Although our observations and those of iBouwman et al.l (120081 ) 
are probing dif f erent regions, it is clear that the two disk layers are connected. Figure 7 of 
Lommen et al.l (120071 ) shows that there is a correlation between the strength/shape of the 
10 /im feature and the slope, a, of the submillimeter emission. Observations of PDS 66 
are consistent with this correlation, with a = —2.4, and the normalized Fu.s/Fg.s ~ 0.97 
and peak -Fio/xm ~ 1.5 (IBouwman et alJl2008h . Although the statistics are still limited, the 
interpretation of this trend is that there is a depletion of sub-/im sized particles because they 
stick together and eventually grow to larger particles. During this grain growth process, the 
10 /im feature flattens out and the slope of the millimeter spectrum becomes shallower, due 
to the changing opacity. We expect that larger grains will s ettle to the midp l ane, c ausing 
the disk to appear less flared. This is exactly what is seen by IBouwman et al.l (120081 ) (their 
Figure 6), where the slope of the FIR emission is correlated with the strength of the 10/iin 
feature (or, as just mentioned, the slope in the submillimeter). Since the slope of the FIR 
emission decreases as the disk becomes less flared, it is clear that increased grain growth, 
traced by either the 10/im feature or the millimeter slope, leads to more settling, and a less 
flared disk. Therefore, o ur observation of a steep er mm-slope is consistent with the flattened 
10 /im feature found by IBouwman et al.l (120081 ). this is most likely accompanied by a less 
flared disk geometry. 
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5.2. Planet-Forming Potential of PDS 66 



It is not clear if giant planets form within all disks, since current statistics are still 
limited. At le ast 80% of the youngest pre-main sequence stars are thought to ha ve circum- 
stellar disks ([Hillenbrand et al.l Il998l ; lHaisch et al.l l200ll ; iHernandez et al.l 120071). whereas 



only 10.5% of stars have gas giant planets with orbital periods <2000 days (jCumming et al. 



20081 ). One prerequisite to form giant planets by core accretion is sufficient disk mass in 
the region where they form, whi ch takes precedence over the shape of the surface density 
profile (e.g. iRaymond et al.ll2005l ). On the theoretical side, core accretion models of planet 
formation require <5 Myr to form Jupiter with a core mass ranging between 5-16 M e (e.g 
Hubickyj et al. II2005I ). If this timescale is right, then long-lived disks such as PDS 66 have 
had enough time to form a gi ant planet. Higher m ass stars have a greater chance of har- 
boring planets at larger radii (1 Johnson et al.l 120071 ) . perhaps simply because more material 
i s available; and su pporting observations suggest that disk mass increases with stellar mass 
(INatta et al.ll2000l ). Could the disk of PDS 66 have enough mass to form planets? 



Pascucci et al.l (120071 ) observed Nell emission towards PDS 66 and that the star is 



actively accreting, providing ample evidence that there is still gas in the disk. We turn 
to our mm data to address how much material may be available to form a planet. As 
Andrews fc Williams I (120071 ) point out, the surface density constrained by observations is 
merely a snapshot in time, while the formation of planets is more likely dependent on the 
evolution of the disk mass (or surface density) in the region where planets might form. 
If the estimated range in the surface density at 5 AU in §4.21 are correct, we can com- 
pare that to_&e minimum models (e.g 



Boss 



1998 



2005 



Pol 



ack et all 11999 : llnaba et al.l 12003 : iDurisen et all 120051 ). Figure 17 of 



Andrews fc Williams I (120071 ) shows the cumulative mass requirements to form giant planets 
from current models. We also find that our upper limit on the total disk mass that is too 
low to form giant planets by about an order of magnitude from the minimum mass require- 
ment. We cannot rule out, however, that a planet may have already formed, though there 
is no evidence of gaps having developed in the SED associated with a forming giant planet. 
A limitation of our analysis is that our unresolved millimeter observations are integrated 
over the disk, and so cannot constrain the mass in potential planet forming regions of the 
disk. Future observations offering significant improvement in spatial resolution in the posited 
planet-forming region, such as with ALMA, will improve constraints on the surface density 
in the regions where planets would form. 
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6. Summary 

We have carried out a detailed analysis of the ~13 Myr old cTTS, PDS 66, using data 
from 0.1 /im out to 12 mm. We conclude the following: 

1. Analysis of the HST/NICMOS scattered light imagery reveals a near-IR surface 
brightness radial distribution going as r -4 53 , detected out to 2" or 170 AU. The inclination 
(32±5°) was derived from the asymmetry in the disk brightness, assuming a mean axial 
ratio=l:0.848. These parameters can be used to remove ambiguity in future modeling efforts 
that could be partially constrained by the surface brightness profile. 

2. The power-law on the opacity spectrum is f3 < 0.5, which was derived from the slope 
of the millimeter spectrum and reflects a maximized correction for optically thick emission. 
The low value of j3 suggests grain growth from the assumed initial ISM sizes, and theoretical 
models suggest that /3 <1 implies millimeter sized grains. 

3. Using a single-temperature blackbody (assuming the inclination from our scattered 
light imaging), we fit the near-IR portion of the SED by varying the inner rim radius and 
temperature. We find that the best-fit parameters are R r i m =0.07 AU and T r j m =1280 K, 
consistent with inner dust disks of other cTTS being truncated at the dust sublimation 
temperature. 

4. Using a geometrically flat, optically thin disk model, we estimate the total surface 
density at 5 AU to 3 < S5 < 11 g/cm 2 , which leads to constraints on the total disk mass 
of 0.0004< Md < O.OOlMffi. This is in accordance with a previous total (dust+gas) disk 



mass estimate ([Carpenter et al.ll2005l ) of 0.005 M (corresponding to 17 M e in dust and 5.3 
M jupiter in gas, assuming a dust-to-gas ratio of 100). These mass constraints indicate that 
this disk is unlikely to form a Jupiter mass planet now, although we found that variations 
in the surface density profile or opacity can lead to disk masses that are higher by as much 
as a factor of 5. This result does not rule out that the possibility that a planet may already 
have formed, though there is no evidence to support this. 
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